Information flow during gene activation by signaling molecules: ethylene transduction in Arabidopsis cells as a study system

<p>Abstract</p> <p>Background</p> <p>We study root cells from the model plant <it>Arabidopsis thaliana </it>and the communication channel conformed by the ethylene signal transduction pathway. A basic equation taken from our previous work relates the probability of expression of the gene <it>ERF</it>1 to the concentration of ethylene.</p> <p>Results</p> <p>The above equation is used to compute the Shannon entropy (<it>H</it>) or degree of uncertainty that the genetic machinery has during the decoding of the message encoded by the ethylene specific receptors embedded in the endoplasmic reticulum membrane and transmitted into the nucleus by the ethylene signaling pathway. We show that the amount of information associated with the expression of the master gene <it>ERF</it>1 (Ethylene Response Factor 1) can be computed. Then we examine the system response to sinusoidal input signals with varying frequencies to determine if the cell can distinguish between different regimes of information flow from the environment. Our results demonstrate that the amount of information managed by the root cell can be correlated with the frequency of the input signal.</p> <p>Conclusion</p> <p>The ethylene signaling pathway cuts off very low and very high frequencies, allowing a window of frequency response in which the nucleus reads the incoming message as a sinusoidal input. Out of this window the nucleus reads the input message as an approximately non-varying one. From this frequency response analysis we estimate: a) the gain of the system during the synthesis of the protein ERF1 (~-5.6 dB); b) the rate of information transfer (0.003 bits) during the transport of each new ERF1 molecule into the nucleus and c) the time of synthesis of each new ERF1 molecule (~21.3 s). Finally, we demonstrate that in the case of the system of a single master gene (<it>ERF</it>1) and a single slave gene (<it>HLS</it>1), the total Shannon entropy is completely determined by the uncertainty associated with the expression of the master gene. A second proposition shows that the Shannon entropy associated with the expression of the <it>HLS</it>1 gene determines the information content of the system that is related to the interaction of the antagonistic genes <it>ARF</it>1, 2 and <it>HLS</it>1.</p

Flower Development as an Interplay between Dynamical Physical Fields and Genetic Networks

In this paper we propose a model to describe the mechanisms by which undifferentiated cells attain gene configurations underlying cell fate determination during morphogenesis. Despite the complicated mechanisms that surely intervene in this process, it is clear that the fundamental fact is that cells obtain spatial and temporal information that bias their destiny. Our main hypothesis assumes that there is at least one macroscopic field that breaks the symmetry of space at a given time. This field provides the information required for the process of cell differentiation to occur by being dynamically coupled to a signal transduction mechanism that, in turn, acts directly upon the gene regulatory network (GRN) underlying cell-fate decisions within cells. We illustrate and test our proposal with a GRN model grounded on experimental data for cell fate specification during organ formation in early Arabidopsis thaliana flower development. We show that our model is able to recover the multigene configurations characteristic of sepal, petal, stamen and carpel primordial cells arranged in concentric rings, in a similar pattern to that observed during actual floral organ determination. Such pattern is robust to alterations of the model parameters and simulated failures predict altered spatio-temporal patterns that mimic those described for several mutants. Furthermore, simulated alterations in the physical fields predict a pattern equivalent to that found in Lacandonia schismatica, the only flowering species with central stamens surrounded by carpels

Continuous-time modeling of cell fate determination in Arabidopsis flowers

<p>Abstract</p> <p>Background</p> <p>The genetic control of floral organ specification is currently being investigated by various approaches, both experimentally and through modeling. Models and simulations have mostly involved boolean or related methods, and so far a quantitative, continuous-time approach has not been explored.</p> <p>Results</p> <p>We propose an ordinary differential equation (ODE) model that describes the gene expression dynamics of a gene regulatory network that controls floral organ formation in the model plant <it>Arabidopsis thaliana</it>. In this model, the dimerization of MADS-box transcription factors is incorporated explicitly. The unknown parameters are estimated from (known) experimental expression data. The model is validated by simulation studies of known mutant plants.</p> <p>Conclusions</p> <p>The proposed model gives realistic predictions with respect to independent mutation data. A simulation study is carried out to predict the effects of a new type of mutation that has so far not been made in <it>Arabidopsis</it>, but that could be used as a severe test of the validity of the model. According to our predictions, the role of dimers is surprisingly important. Moreover, the functional loss of any dimer leads to one or more phenotypic alterations.</p

Vertebrate Paralogous MEF2 Genes: Origin, Conservation, and Evolution

BACKGROUND: The myocyte enhancer factor 2 (MEF2) gene family is broadly expressed during the development and maintenance of muscle cells. Although a great deal has been elucidated concerning MEF2 transcription factors' regulation of specific gene expression in diverse programs and adaptive responses, little is known about the origin and evolution of the four members of the MEF2 gene family in vertebrates. METHODOLOGY/PRINCIPAL FINDINGS: By phylogenetic analyses, we investigated the origin, conservation, and evolution of the four MEF2 genes. First, among the four MEF2 paralogous branches, MEF2B is clearly distant from the other three branches in vertebrates, mainly because it lacks the HJURP_C (Holliday junction recognition protein C-terminal) region. Second, three duplication events might have occurred to produce the four MEF2 paralogous genes and the latest duplication event occurred near the origin of vertebrates producing MEF2A and MEF2C. Third, the ratio (K(a)/K(s)) of non-synonymous to synonymous nucleotide substitution rates showed that MEF2B evolves faster than the other three MEF2 proteins despite purifying selection on all of the four MEF2 branches. Moreover, a pair model of M0 versus M3 showed that variable selection exists among MEF2 proteins, and branch-site analysis presented that sites 53 and 64 along the MEF2B branch are under positive selection. Finally, and interestingly, substitution rates showed that type II MADS genes (i.e., MEF2-like genes) evolve as slowly as type I MADS genes (i.e., SRF-like genes) in animals, which is inconsistent with the fact that type II MADS genes evolve much slower than type I MADS genes in plants. CONCLUSION: Our findings shed light on the relationship of MEF2A, B, C, and D with functional conservation and evolution in vertebrates. This study provides a rationale for future experimental design to investigate distinct but overlapping regulatory roles of the four MEF2 genes in various tissues

Evolving Sensitivity Balances Boolean Networks

We investigate the sensitivity of Boolean Networks (BNs) to mutations. We are interested in Boolean Networks as a model of Gene Regulatory Networks (GRNs). We adopt Ribeiro and Kauffman’s Ergodic Set and use it to study the long term dynamics of a BN. We define the sensitivity of a BN to be the mean change in its Ergodic Set structure under all possible loss of interaction mutations. Insilico experiments were used to selectively evolve BNs for sensitivity to losing interactions. We find that maximum sensitivity was often achievable and resulted in the BNs becoming topologically balanced, i.e. they evolve towards network structures in which they have a similar number of inhibitory and excitatory interactions. In terms of the dynamics, the dominant sensitivity strategy that evolved was to build BNs with Ergodic Sets dominated by a single long limit cycle which is easily destabilised by mutations. We discuss the relevance of our findings in the context of Stem Cell Differentiation and propose a relationship between pluripotent stem cells and our evolved sensitive networks

Identification and Characterization of a Mef2 Transcriptional Activator in Schistosome Parasites

Myocyte enhancer factor 2 protein (Mef2) is an evolutionarily conserved activator of transcription that is critical to induce and control complex processes in myogenesis and neurogenesis in vertebrates and insects, and osteogenesis in vertebrates. In Drosophila, Mef2 null mutants are unable to produce differentiated muscle cells, and in vertebrates, Mef2 mutants are embryonic lethal. Schistosome worms are responsible for over 200 million cases of schistosomiasis globally, but little is known about early development of schistosome parasites after infecting a vertebrate host. Understanding basic schistosome development could be crucial to delineating potential drug targets. Here, we identify and characterize Mef2 from the schistosome worm Schistosoma mansoni (SmMef2). We initially identified SmMef2 as a homolog to the yeast Mef2 homolog, Resistance to Lethality of MKK1P386 overexpression (Rlm1), and we show that SmMef2 is homologous to conserved Mef2 family proteins. Using a genetics approach, we demonstrate that SmMef2 is a transactivator that can induce transcription of four separate heterologous reporter genes by yeast one-hybrid analysis. We also show that Mef2 is expressed during several stages of schistosome development by quantitative PCR and that it can bind to conserved Mef2 DNA consensus binding sequences

Identification and Characterization of RcMADS1, an AGL24 Ortholog from the Holoparasitic Plant Rafflesia cantleyi Solms-Laubach (Rafflesiaceae)

10.1371/journal.pone.0067243PLoS ONE86-POLN

Noise and Robustness in Phyllotaxis

A striking feature of vascular plants is the regular arrangement of lateral organs on the stem, known as phyllotaxis. The most common phyllotactic patterns can be described using spirals, numbers from the Fibonacci sequence and the golden angle. This rich mathematical structure, along with the experimental reproduction of phyllotactic spirals in physical systems, has led to a view of phyllotaxis focusing on regularity. However all organisms are affected by natural stochastic variability, raising questions about the effect of this variability on phyllotaxis and the achievement of such regular patterns. Here we address these questions theoretically using a dynamical system of interacting sources of inhibitory field. Previous work has shown that phyllotaxis can emerge deterministically from the self-organization of such sources and that inhibition is primarily mediated by the depletion of the plant hormone auxin through polarized transport. We incorporated stochasticity in the model and found three main classes of defects in spiral phyllotaxis – the reversal of the handedness of spirals, the concomitant initiation of organs and the occurrence of distichous angles – and we investigated whether a secondary inhibitory field filters out defects. Our results are consistent with available experimental data and yield a prediction of the main source of stochasticity during organogenesis. Our model can be related to cellular parameters and thus provides a framework for the analysis of phyllotactic mutants at both cellular and tissular levels. We propose that secondary fields associated with organogenesis, such as other biochemical signals or mechanical forces, are important for the robustness of phyllotaxis. More generally, our work sheds light on how a target pattern can be achieved within a noisy background

Comparison of evolutionary algorithms in gene regulatory network model inference

Background: The evolution of high throughput technologies that measure gene expression levels has created a data base for inferring GRNs (a process also known as reverse engineering of GRNs). However, the nature of these data has made this process very di±cult. At the moment, several methods of discovering qualitative causal relationships between genes with high accuracy from microarray data exist, but large scale quantitative analysis on real biological datasets cannot be performed, to date, as existing approaches are not suitable for real microarray data which are noisy and insu±cient. Results: This paper performs an analysis of several existing evolutionary algorithms for quantitative gene regulatory network modelling. The aim is to present the techniques used and o®er a comprehensive comparison of approaches, under a common framework. Algorithms are applied to both synthetic and real gene expression data from DNA microarrays, and ability to reproduce biological behaviour, scalability and robustness to noise are assessed and compared. Conclusions: Presented is a comparison framework for assessment of evolutionary algorithms, used to infer gene regulatory networks. Promising methods are identi¯ed and a platform for development of appropriate model formalisms is established